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Abstract 

We study in detail coherent phase control of femtosecond resonance-mediated (2+1) three-photon 
absorption and its dependence on the spectral bandwidth of the excitation pulse. The regime is the 
weak- field regime of third perturbative order. The corresponding interference mechanism involves a 
group of three-photon excitation pathways that are on resonance with the intermediate state and a 
group of three-photon excitation pathways that are near resonant with it. The model system of the 
study is atomic sodium (Na), for which experimental and numerical-theoretical results are obtained. 
Prominent among the results is our finding that with simple proper pulse shaping an increase in 
the excitation bandwidth leads to a corresponding increase in the enhancement of the three-photon 
absorption over the absorption induced by the (unshaped) transform-limited pulse. For example, 
here, a 40-nm bandwidth leads to an order-of-magnitude enhancement over the transform-limited 
absorption. 
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I. INTRODUCTION 



Coherent control utilizes the coherent nature of the light in order to control transition 
probabilities of the irradiated quantum system to desired target states [H El El IH E] ■ Due 
to their broad spectrum, femtosecond pulses provide a wide coherent band of such photo- 
induced pathways. The control is achieved by manipulating interferences among multiple 
initial-to-final state-to-state pathways. Constructive interferences enhance the transition 
probability, while destructive interferences attenuate it. The interferences are manipulated 
by shaping the pulse, i.e., manipulating the phase, amplitude, and/or polarization of its 
various spectral components. 

Prominent among the processes, over which such femtosecond control has been demon- 
strated to be very effective, are multiphoton processes that are of fundamental as well as 
applicative importance for nonlinear spectroscopy and microscopy O El [71 [H |9l [TOl [TTl [121 
[T3l [HI [151 HSl [HI mi nil [201 [21] . For a given energy, manipulating the pulse shape manip- 
ulates the cross-section of the corresponding process. For rational shaping of femtosecond 
pulses, which are favorable for a desired result, the most convenient excitation regime is the 
weak-field regime where the photo-excitation is validly described within the time- dependent 
perturbation theory of the lowest non-vanishing order. For N-photon absorption processes it 
is the N*'^ order. Such a description allows a transformation of the photo-excitation picture 
into the frequency domain, where the interfering pathways and their interference mechanism 
can be identified and serve as a basis for the rational pulse shaping. Corresponding examples 
of past control studies include the control over non-resonant two-photon absorption P, [7j 
and resonance-mediated (1+1) two-photon absorption [8l[9l[Tn] in atoms, over non-resonant 
two-photon absorption and non- resonant three-photon absorption in large molecules [TT] , 
over molecular Raman transitions [I2l, and over molecular coherent anti-Stokes Raman 
scattering (CARS) ^EM- 

In addition to these multiphoton control works, in a recent study we firstly analyzed 
theoretically and demonstrated experimentally coherent phase control over the process of 
atomic resonance-mediated (2+1) three-photon absorption [15]. The present work is a di- 
rect continuation of this first work. It extends the analysis of the corresponding interference 
mechanism and obtains new physical insight regarding its dependence on the spectral band- 
width of the excitation pulse. The model system of the study is atomic sodium (Na). 
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Section |TT] of the paper gives the corresponding frequency-domain theoretical description of 
the process as obtained based on 3''''^-order time-dependent perturbation theory. The paper 



ends with conclusions in Sec.M Section III A describes the Na model system and the families 



of shaped pulses that serve as test cases. It follows by Sec. |III B| that presents corresponding 
experimental results that re-confirm the theoretical description and corresponding numeri- 
cal calculations. Then, Sec. [Ill C presents extended numerical-theoretical results that fully 



reveal the spectral-bandwidth dependence of the process, as discussed in Sec. IV 



II. FREQUENCY-DOMAIN THEORETICAL DESCRIPTION 

As shown in Fig. [T| the resonance-mediated (2+1) three-photon absorption process con- 
sidered here involves an initial state \g), a final state |/), and an intermediate state |r). 
The \g) and |r) states are coupled by a non-resonant two-photon coupling provided by a 
manifold of states \v) that are far from resonance, i.e., half the two-photon transition fre- 
quency iOrgjl falls within the pulse spectrum while the \g)-\v) transition frequencies lo^g fall 
outside the pulse spectrum. The coupling is a resonant one-photon coupling, i.e., the 

corresponding transition frequency LOj^. falls within the pulse spectrum. 

In a weak-field regime the three-photon absorption is validly described by third-order 
time- dependent perturbation theory. Accordingly, the time-dependent (complex) amplitude 
a/(t) of state |/) at time t, due to the interaction with the shaped temporal electric field 
e{t), is given by 

a/(t) = -^^f^frfJ'rvfJ'vg [ [ / e{ti)e{t2)e{t3) 

X exp{iLJfrti) exp{iLJrvt2) exp{iuygt3)dtidt2dt2, (1) 

where fimn = (^| A* \n) and Umn = {E^ — En)/h are, respectively, the dipole matrix element 
and the transition frequency between a pair of states |m) and \n). The transition frequencies 
ojrg, ujfr, and Ufg satisfy the relation Ufg = Ufj. + uj^g. 

Within the frequency-domain framework, the spectral field of the pulse E{uj) = 
\E{uj)\exp[i^{uj)] is given as the Fourier transform of e{t), with |-E'(co')| and ^{u) being, 
respectively, the spectral amplitude and phase of frequency u. For the unshaped transform- 
limited (TL) pulse, $(co') = for any u. We also define the normalized spectral field 
E{ij) = E{ijj)/ \Eq\ = E{uj) exp [^$(ci;)] that represents the pulse shape, where \Eq\ is the 
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peak spectral amplitude. This allows to clearly distinguish in the expressions given below 
between the dependence on the pulse intensity and the dependence on the pulse shape. 
The maximal spectral intensity Jq is proportional to l-EoP (-^o oc |-EoP)- The resonance- 
mediated (2+1) three-photon excitation scheme implies that E{uJrg/2) 



7^0, 



0, 



and 



E{ujfr) 7^ 0. Then, the final amplitude Aj'^^^ = cif{t — > oo) of state |/) after the pulse 



is over {t oo) is obtained to be [15] 

.(2+1) _ 1 ,,2 17-1 |3 \ A{2+l)on-res . (2+l)near-resl /^n 

^/ - ^/^/r/^rs Rol + Aj J, {2) 

where fifr is the dipole matrix element and fi'^g is the \r)-\g) effective non-resonant 

two-photon dipole coupling The terms '^^'^ and '"'^^ qj-q given by 

^(2+l)on-res ^ E {cU fr) A^'\u;rg) , (3) 

/oo 1 ^ 
|A(2)K,-5)i?(^;. + 5)rf5, (4) 

with 

/oo ^ ^ 
E{iu)E{n - uj)duj . (5) 
-oo 

The expressions of Eqs. (|2])-([5]) are equivalent to the ones given in Ref. [I5]. However, as 
shown below, they are given here in a different form that provides deeper and more intuitive 
physical insight to the femtosecond excitation process. 

The final amplitude of |/) coherently interferes all the possible three-photon path- 

ways from \g) to |/), i.e., coherently integrates over all their corresponding amplitudes. 
Several pathway examples are shown schematically in Fig. [T} Each three-photon pathway is 
either on-resonance or near-resonance with the intermediate state |r), with a corresponding 
detuning 6. A resonance-mediated (2+1) three-photon pathway corresponding to a detuning 
6 involves a non-resonant absorption of two photons to a total excitation energy of uJrg — ^ 
and the absorption of a third photon of frequency luj^ + 6. The term j^^gj^_ 
feres all the on-resonant pathways (5 = 0), while the term ^(2+i)"ear--res j^^^gj^fgj^gg g^jj ^j-^g 
near-resonant pathways {6 ^ 0) with an amplitude weighting of 1/5 (i.e., inversely propor- 
tional to its detuning). The on-resonant pathways are excluded from ^(2+i)near res 
Cauchy's principal value operator p. Thus, the resulting absorption is determined by both 
intra-group and inter-group interferences involving these two groups of three-photon exci- 
tation pathways. Both g^j^^ j^(2+i)near-res ^^^.j^^ g^j,g exprcsscd using the term 
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A^"^") (il) that is proportional to the non-resonant two-photon absorption amphtude to a total 
transition frequency of Q. It interferes all the corresponding two-photon pathways. The 
phase associated with each two-photon pathway is equal to the sum of the phases of the 
two corresponding photons. So, with the transform-limited (TL) pulse all the two-photon 
pathways acquire zero phase for any value of Q and, thus, interfere one with the other in a 
fully constructive way. For a given field spectrum E{uj) , this leads to the maximal A^'^\il,) 
for any value of Q. For a given Q, any spectral phase pattern that is anti-symmetric around 
Q/2 also induces fully constructive interferences among all the two-photon pathways of Q 
transition frequency, and thus it induces the corresponding maximal A'^'^\^l) as the TL pulse. 
As can be seen from Eq. (pj), A^'^\n) is actually the (complex) convolution function of the 



(complex) field E{uj) evaluated at fi. Since y^(2+i)on-res proportional to A^'^\ujrg), it is 
actually proportional to the weak-field non-resonant two-photon absorption amplitude to 
|r), as obtained by second-order perturbation theory |6j. 



.(2+1) 2 



serves as the measure 



The final population of state |/), given by Pj- 

for the resonance-mediated (2-1-1) three-photon absorption. All the results presented here 
are given relative to the absorption induced by the transform-limited (TL) pulse. So, for 
convenience, we introduce the TL-normalized three-photon absorption measure pj^~^^^ = 
p{2+i) I p{2+i) ^ ^^^^^ -^/TL^ i^ final I/) population induced by the TL pulse. 



III. RESULTS 



A. Model system and pulse-shape test cases 

The model system of the study is the sodium (Na) atom |23], with the corresponding 
excitation scheme shown in Fig. [l] It involves the 3s ground state as \g)^ the 4s state as |r), 
and the 7p state as |/). The transition frequency u^g = uj4^s,3s = 25740 cm^-*^ corresponds to 
two 777 nm photons and the transition frequency ujfr = uj-ip^^g = 12801 cm~^ corresponds 
to a one 781.2 nm photon. The 3s-4s non-resonant two-photon coupling originates from the 
manifold of p-states, particularly from the 3p state [uj^p^^s ~ 16978 cm~^ (589 nm)]. The 
absorption measure is the final 7p population pj^^^-* = Py^^^^ 

As the test cases for coherent phase control, the present study uses shaped femtosecond 
pulses having two sets of spectral phase patterns, <|)(^'^'^*''p)(co') and <|)(^^'^*''p)(ci;). They are 
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shown schematically in the inset of Fig. [l| The <|)(s^'^**^p)(u;) set includes the phase patterns 
of a single vr step, each characterized by the vr step position oo^^^p 

includes the phase patterns composed of two vr steps, with the two steps located at equal 
distance along the frequency axis from both sides (i.e., to the blue and to the red) of half 
the two-photon transition frequency u;4s,3s/2, i.e., ^^4.,3s/2 -t^St^s^p^ = ^Sht-S -^4.,3s/2 = 



Ag^p***^''\ In our previous work [13] the single-step phase patterns have been shown to be 
very effective in controlling the resonance- mediated (2+1) three-photon absorption. The 
double-step phase patterns are additionally included here since, as explained below, they all 

1(24 



induce the same maximal on-resonant component A^^^^^"" while each of them induces 



a different near-resonant component ^(2+i)near res 



^7p 



The vr amplitude of all the phase steps, i.e., generally ^{uj) = or vr, implies 
that the field E{uj) is always a real (positive or negative) quantity for any u. Thus, 
the corresponding Aj^ 



(2+l)on—res ^ t (2+l)near—res i_- i i • • i 

and Ar^„ are, respectively, purely imagmary and 



real [see Eqs. (3)-(4)], and P^^ 

2 



^7p 



(2+1) 



(2+1) 



7p 

and 



A 



{2+l)on~res 
7p 



+ 



A 



A 



'{2+l)near—res 
7p 



A 



oc 

(2+l)near—res 
2 



A 



{2+l)on—res 
7p 



+ 



A 



{2+l)near—res 
7p 



7p 



(2+l)on— res 
7p 



A 



(2+l)o7i-re, 
7p 



Consequently, 
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IP, 



(2+1) 
7p,TL 



where A 

(2+i)near res j pjj^'t^^. In other words, with the spectral phase pat- 
terns discussed here, the absorption is determined only by intra-group interferences taking 
place separately within each of the on-resonant and near-resonant excitation groups, without 
any inter-group interferences taking place. This fact greatly simplify, without reduction of 
generality, the ability to gain physical insight into the pulse-bandwidth dependence of the 
resonance-mediated (2+1) three-photon absorption process. 



B. Experimental vs. numerical-theoretical results 

This section presents experimental results for the Na atom, re-confirming the above the- 
oretical description with spectral parameters, i.e., spectral bandwidth and the above phase 
patterns, that are different from those used in our previous work [15] . Experimentally, atomic 
sodium vapor is produced in a static chamber at 300°C (Na partial pressure of ~0.1 Torr) 
with 10- Torr Ar buffer gas. It is irradiated at a 1-kHz repetition rate with linearly-polarized 
shaped femtosecond laser pulses of 9.8-nm (160-cm~^) intensity bandwidth (FWHM) around 
781 nm (12804 cm~^). The spectrum shape is a modified Gaussian with a slight asymmetry 
toward low frequencies. It is shown in the inset of Fig. [2Fal). Experiments are also con- 
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ducted with a spectrum that results from blocking about half of this full spectrum at a cutoff 
wavelength of 783 nm (12771 cm^^). It is shown in the inset of Fig. [2|^bl). In both spec- 
tral cases, the temporal peak intensity of the corresponding transform-limited (TL) pulse 
is below 10^ W/cm^. The laser pulses undergo shaping in a 4/ optical setup incorporating 
a pixelated liquid-crystal spatial light phase modulator [22]. The effective spectral shaping 
resolution is Sughaping = 2.05 cm~^ (0.125 nm) per pixel. Following the interaction with a 
pulse, the Na population excited to the 7p state undergoes radiative and collisional decay 
to lower excited states, including the Ad, 5d, 6d, and 6s states. The fluorescence emitted 
in their decay to the 3p state serves as the relative measure for the total final 7p popu- 
lation Pjp^^\ It is optically measured at 90° to the beam propagation direction using a 
spectrometer coupled to a time-gated camera system. 

Figure [2] compares the experimental results (circles or squares) with corresponding theo- 
retical results (solid lines) for the resonance- mediated (2+1) three- photon absorption in Na. 
The theoretical results are calculated numerically using Eqs. ([2])-(|5|, using a grid with a bin 
size equal to the experimental spectral shaping resolution 6u shaping- Figures |2{al) and|2|(a2) 
(first-row panels) correspond to the full pulse spectrum [inset of Fig.[2]^al)], while Figs. |2]^bl) 
and|2](b2) (second-row panels) correspond to the blocked pulse spectrum [inset of Fig.|2](bl)]. 
Figures [2]^al) and[2]^bl) (left-column panels) present results for the single-step set of spectral 
phase patterns <|)(^^'^*'=p)(co'). The corresponding traces show the TL-normalized final 7p pop- 
ulation -Pyp'*'^'' (see above) as a function of the single vr-step position cUgf^p^*'^'''' . Figures 2|a2) 
and|2](b2) (right-column panels) present results for the double-step set of spectral phase pat- 
terns <|)(°^'^*'^p)(co'). The corresponding traces show the TL-normalized final 7p population 
Pjp^^^ (see above) as a function of the left vr-step position c^j^^step^- 

As can be seen, there is an excellent quantitative agreement between the experimental 
results ("real experiment") and the numerical-theoretical results ("computer experiment"). 
Hence, supported also by our previous work [12] , the validity and accuracy of the theoretical 
description and numerical calculations are re-confirmed. As such, in the following they 
serve us for the detailed study of the vr-traces, their spectral-bandwidth dependence, and 
the corresponding interference mechanism. 
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C. Numerical-theoretical results 



Figure [3] presents theoretical results (solid lines) that have been calculated numerically 
for the resonance-mediated (2+1) three-photon absorption in Na with excitation pulse of 
different spectral bandwidth. In all the cases the intensity spectrum is a Gaussian centered 
at 780 nm (12821 cm~^). The corresponding bandwidth (FWHM) values (Acu) are (a) 
5 nm (82 cm"-"^), (b) 9 nm (148 cm^-"^), (c) 15 nm (247 cm^-*^), (d) 25 nm (411 cm"^), and 
(e) 40 nm (658 cm~^). Each row in the figure corresponds to a different bandwidth. As 
indicated above, the numerical calculations are based on Eqs. (|2])-([5]) and use a grid with a 
bin size equal to the experimental spectral shaping resolution. 

The left-column panels of Fig. [3] correspond to the single-step spectral phase patterns 
^(SNstep)^^^_ In solid lines they show the TL-normalized final 7p population -Pyp"*^^^ as a 
function of the single vr-step position co'gfgp^*'^^^ Additionally, they also show the correspond- 

(dotted lines) and the near-resonant component 



ing on-resonant component 

2 

j^{2+l)near—res 



A (2+l)ora— res 
^7p 



A (2+l)on— res 
^7p 



+ 



A {2+l)near—res 
^7p 



see 



(dashed lines), satisfying P^^^^^ 
above). As can be seen, the shape of the traces strongly depends on the spectral band- 
width. For any given bandwidth, the three-photon absorption is tunable over a wide range 
of values, from a low attenuated level below 0.01 of the TL absorption to a high enhanced 
level with a value that strongly depends on the spectral bandwidth. Whenever this en- 
hanced absorption exceeds the TL absorption, it occurs when cUsfep'^*"'^'' = ^ipM^ i-^-' when 
the vr-step is located at the 7p-4s transition frequency. This is actually the location of the 
prominent feature of the single-step traces. As seen in Fig. [3} as the spectral bandwidth 



increases from 5 nm to 9 nm to 15 nm, the value of -P^p'^^'* at ti^step'''''^'' = ^ipAs decreases, 
respectively, from 2.5 to 1.1 to below 0.01. Then, further increase in the bandwidth leads to 
its increase to a value of 1.4 at 25-nm bandwidth and to a value of 4.1 at 40-nm bandwidth. 
Another prominent bandwidth dependence is observed for the TL pulse regarding the rela- 
tive weight of the corresponding on-resonant and near-resonant amplitude components, as 
it is reflected in the asymptotes of the single-step traces: As the bandwidth increases, the 
TL on-resonant component dominates more and more the total absorption, approaching the 

2 



(SNstep) 



limit of 



A 



{2+l)near—Tes 
7p,TL 



A 



(2+l)on— res 
7p,TL 



and 



A 



{2+l)on— res 
7p,TL 



The right-column panels of Fig. [3] correspond to the double-step spectral phase patterns 

.(2- 

7p 



^(DBstep)^^^_ They show the TL-normalized final 7p population P^t^^^ as a function of the left 
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TT-step position ^^ic'ft^step "* • C)ne should note that the maximal value of the y-axis in these traces 
is higher than in the traces of the single-step phase patterns (the left-column panels). The 
prominent feature of the double-step traces is the strong absorption enhancement occurring 
when i^le^t^step^ ~ ^ipAs- This is also the step-position region, where the dependence on 
the spectral bandwidth is the strongest. The corresponding enhancement value is always 
above one, i.e., above the TL absorption, and it continuously increases as the bandwidth 
increases. It starts from a value of 2.8 at 5-nm bandwidth [Fig. |3]^a2)] and reaches a value 
of 7.2 at 40- nm bandwidth [Fig. |3|e2)]. When i^ie^t^stTp'' is far from t<;7p,4s, the three-photon 
absorption is always kept on a level that is very close to or equal to the TL absorption. The 
only exception is the case of 5-nm bandwidth [Fig. |3]^a2)], where at step positions higher 
than uj7p^4s the absorption is reduced to a level of about half the TL absorption. 

As can be seen in Fig. 3 , the main features of the traces when cUgfep'^*'^^^ = t^7p,4s and when 

(DBstep) • • J. r J.U j. j. < (2+l)near-res rpu j. 

"^left-step ~ ^TpAs origmatc from the near-resonant component Aj^ . ihus, to 

complete the picture, Fig.|4]presents extended results for the spectral-bandwidth dependence 
of the resonance-mediated (2-1-1) three-photon absorption when cUgfgp'^*'''''' = oajp^^s in a single- 
step pattern [Fig. 4[a)] and when (^[^f^.^l^p^ = ujjp^^s in a double-step pattern [Fig. 4[b)]. 

i_r ' 2 '-'^ 



The figure presents the corresponding values of pl^"*"^^ (solid hues), A\V''''°'^ (dotted 



/ 



j'(2+l)on— res 



lines), and A^p"'"^^"'^"'^ (dashed lines) as a function of the spectral bandwidth. The 
bandwidth values span the wide range of 2 to 40 nm that corresponds, respectively, to a TL 
pulse duration of 450 down to 22 fs. As can be seen, with the single-step phase pattern of 
djgfep'^*'^^'' = ojjp^As, the degree of absorption is between 0.01 to 4 times the TL absorption, i.e., 
it is either enhanced or attenuated as compared to the TL case. There is a local peak of value 
2.55 at a bandwidth of 5.25 nm, and a minimum of zero theoretical value (i.e., three-photon 
dark pulse) at a bandwidth value of 14.6 nm. Then, above 14.6 nm, pj^^^-* is monotonically 
increases with the increase in the bandwidth. At 23-nm bandwidth it reaches a value of 1, 
at 32 nm it reaches again a value of 2.5, and at 40 nm it reaches a value of 4.1 . On the 
other hand, with the double-step phase pattern of c^ieft^step'' ~ ^7p,4s; the absorption is always 
above the TL absorption and the corresponding enhancement monotonically increases with 
the increase in the spectral bandwidth. At 2-nm bandwidth Pyp"*"^^ is of a value of 1, and it 
reaches a value of 7.2 at a bandwidth of 40 nm. As can be seen in Fig. |4| in both these single- 
step and double-step patterns having the step at oj-jp^/^s, over the full bandwidth range, the 
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near-resonant component 

2 



A 



{2+l)near—res 
7p 



is much larger than the on-resonant component 



A 



{2+l)on—res 
7p 



IV. ANALYSIS AND DISCUSSION 



Below, the main features of the results presented above are analyzed based on the the- 
oretical description formulated in Eqs. ([2])-(|5]). As noted above, for the spectral phase pat- 
terns used here, the resonance-mediated (2+1) three-photon absorption is determined only 
by intra-group interferences within each of the on-resonant and near-resonant excitation 
groups that correspond, respectively, to A'^^^^^""' and A^^^^''^^'^'^ , with the relations 

p(2+l) 



OC 



A {2+l)on—res 
^7p 



+ 



A {2+l)near—res 
^7p 



and 



A (2+l)ora— res 
^7p 



+ 



A {2+l)near—res 
^7p 



The analysis uses the numerical results presented in Fig. |5j Each panel in the figure 
corresponds to a pulse having a different spectral bandwidth for its Gaussian intensity 
spectrum around 780 nm. The corresponding bandwidth values are (a) 5 nm, (b) 9 nm, 
(c) 15 nm, (d) 25 nm, and (e) 40 nm, as they are in the results shown in Fig. [s] Each panel 
shows results for the transform-limited (TL) pulse (dashed black line), for the shaped pulse 
with the single-step phase pattern of cUgtep'^*'''''' = ^7p,4s (solid gray line), and for the shaped 
pulse with the double-step phase pattern of i^l^t^step'' ~ ^7p,4s (solid black line). For the 
different cases, the values of A^'^\il = Uis^ss — [see Eq. ^] are shown as a function of the 
detuning 6 from the 4s state. The detuning (x-axis) values are actually given as normalized 
values of 6/Auj, with Auj being the bandwidth of the intensity spectrum. 

As seen from Eq. (3), the value of the on-resonant component Ayp"*"^^"" ^'^^ is proportional 
to A^'^^^\0) = A^'^\ij4s,3s)E{ujYp^4s) that corresponds to zero detuning 5 = 0. On the other 
hand, as seen from Eq. Q, due to the 1/5- weighting, the result of the p-integration of 
y^{2+i)near-res (dominated by the integration over small non-zero values of |5| around 6=0. 
Due to the sign change of 1/5 for ± |5|, the p- integration result is highly sensitive to the 
symmetry of the integrand term A^'^'^^\6) = A^'^\uj4,s,3s — S)E{uj7p^4s + S) around 6=0, i.e., its 
relative magnitude and relative sign for positive and negative detunings of equal magnitude 
(i.e., ± \6\). The zone of small \6\ around 6=0 is indicated schematically as dashed area in 
Fig.m 
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A. Transform-limited pulse 



First we analyze the case of the transform-hmited (TL) pulse (dashed thin black lines 
in Fig. Since A^'^\^l) is the convolution function of E{uj) evaluated at fi, for the TL 
pulse [^{uj) = 0] it is actually a convolution of E{oo) with a real positive values for any Q. 

(2) 

Thus, with a Gaussian spectrum around uq, A>^l{fl) has a Gaussian shape that is peaked 
at Qpeak = 2ti;o- The corresponding detuning is 6peak = i^4s,3s — Stuo- The width (FWHM) 
of the A!^l(il) Gaussian profile is 2Auj, where Au is the width of the intensity spectrum. 
Since the detuning values given in Fig. [5] are normalized by Au, all the presented TL profiles 
of Ai^l{^l = uJ4s,3s ~ ^) have the same width (of value 2) for any spectral bandwidth Aoo. 
On the other hand, an increase in Au leads to a decrease in the value of Speak/ Au. Hence, 
as seen in Fig. [sj as Auj increases, the peak and the full profile of Ai^l{uj4s,3s — S) shifts to 
smaller values of 6/Auj and the corresponding value of A^^l('^4s,3s)) which determines the 



on-resonant component A 



{2+l)on—res 
7p,TL ' 

I (2) 



mcreases. 



For the TL pulse, both A\,l^{ujis,-is — S) and -Etl (f^Tp^s + S) are real positive quantities for 
any 6. Thus, since the sign of 1/5 is different for and — the amplitudes ^A^/^^^\5) 
contributed by ± \5\ to the p-integration of A^^/^]^^"'^ are of different signs, i.e., they 
interfere destructively one with the other. 

As a result from all the above, an increase of Auj leads to a decrease in the value of 



i 4(2+1)/ 



/ A^/^^^\q) for a given small \5\. Consequently, there is a 
decrease in the relative magnitude between the TL near-resonant component 
and the TL on-resonant component ^7p'tl°" ■ Eventually, at large enough Au;, it leads to 



A 



{2+l)near—res 
7p,TL 



A. 



(2+l)on— res 
7p,TL 



p(2+l) 
' -'^7p,TL 



A {2+l)on—res 
^7p,TL 



and 



i {2+l)on—res 
^7p,TL 



1. As noted 



above, this is indeed the behavior seen in the asymptotes of the single-step traces shown in 
the left-column panels of Fig. [s] (dashed and dotted lines). In other words, for the TL pulse, 
with large enough bandwidth, the absorption amplitudes contributed by positively-detuned 
and negatively-detuned near-resonant three-photon pathways completely cancel out each 
other. Then, the total TL absorption is determined only by the three-photon pathways that 
are on-resonant with the intermediate state. 
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B. Shaped pulse with single-step phase patten of '^gfj^^*''''^ = ^7p,4s 

Next we consider the case of the shaped pulse with the phase pattern of a single vr step at 
SNstcp) ^ ^^^^^^ ^^^^.^ g^^y ^.^^g p.g^ 5). As the corresponding field -EsNstep@a;7p,4s (^^7p,4s + 
6) is a real positive or negative quantity for any u, the corresponding ^sNstepOojyp 4^ ("^^^'Ss ~ ^) 
is also a real quantity that, for a given spectrum, its positive or negative sign depends on 6. 
This is seen in Fig. [5} where the difference from the TL case is clearly observable. Since the 
TL pulse induces fully constructive interferences among all the two-photon pathways that 
contribute to A^'^^Q) for any given Q, with a given spectrum, the TL value of A'^1{uj4s,3s — S) 
is actually the maximal one for any given 6. As such, it serves as an upper limit to the value 

(^) I 

of ^SNstep@a;7p,4s(^4s,3s - 

As can be seen from Fig. jsj the values of ^sNstep®w7p 4^ ('^4s,3s — S) at 6 = and at 5 ~ 
strongly depend on the spectral bandwidth: with 5-nm and 9-nm bandwidth the values for 
5 = and 5 ~ are all positive [panels (a) and (b)], with 15- nm bandwidth the value for 
5 = is almost zero while the values for small ± \6\ are of different signs [panel (c)], and 
with 25-nm and 40-nm bandwidth the values for 5 = and 5^0 are all negative [panels (d) 
and (e)]. On the other hand, with any spectral bandwidth, the field -EsNstep@(^7p,4s (i^7p,4s + 5) 
has different sign for +\6\ and which is in full correlation with the corresponding sign 
of 1/6. Thus, the relative sign between the amplitudes |^sNstep(aa;7p4s('^) contributed by 
± |5| to the p-integration of Afp^sN^Zi^Zis determined only by ^sNstep@a;7p,4s (^4s,3. - 
As such, it depends here on the spectral bandwidth. In Fig. jsj the ^sNstep@LJ7p 4^ ('^4s,3s ~ ^) 
values for ± \6\ are of the same sign, i.e., they interfere constructively one with the other, 
in all the bandwidth cases except for the 15-nm bandwidth, where they are of different 
sign and thus interfere destructively one with the other. The constructive nature of these 
interferences generally leads to the dominance of the corresponding near-resonant component 

^ipS"tep@"J7p.4. on-resonant component Afj^snZprn^rpAs bandwidth values, 

except for values around 15 nm. There, both the near-resonant and on-resonant components 
are very small and have a value close to zero. This is indeed the behavior of the near- resonant 
vs. on-resonant components shown in Fig. |4]^a). 

For analyzing the bandwidth dependence of the absorption relative to the TL pulse, i.e., 
the graph of -Pyp sNstep@a;7p 4^ shown in Fig. 4j(a) (solid line), one also needs to account 

for the case of the TL pulse analyzed in the previous section. At a very narrow bandwidth 
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the phase shaping is expected to have very small or no effect at all on the absorption relative 
to the (unshaped) TL pulse, as is indeed reffected in the value -Prp sNstepOwyp 4^ ~ ^"-"^ 2-nm 
bandwidth. Then, as the bandwidth increases the pulse shaping starts to have an effect. 



For the 5-nm bandwidth, where there is a local peak of -Py^ sNstep' 



that: (i) 
[see Fig. 



< (2+l)on— res 
^7p,SNstep@cJ7p,4, 

^a)], and (ii 



(2+l)on— res 



^7p,TL 

A {2+l)near— res 
^7p,SNstep@i^7p,4s 



St^7p,4s 

since ^SNstep@.;Tp,4.(^4s,3s - 5) 



> 



A 



(2) 
SN; 

(2+l)near— res 
7p,TL 



2.5, one obtains 

(2) 

i A^£(u;4s,3s - 5) 



due to the different nature of the 



interferences between positively- and negatively-detuned excitation pathways for the two 



pulses. Hence, P^pjijltep. 



. (2 + l)on — res 



. (2 + l)near — res 
Tp.SNstopOuyp 4g 



S'^7p,4s 



.(2 + l)on — res 
7p,TL 



. (2 + l)near — res 
7p,TL 



> 1 for 5-nm bandwidth. 



Then, as the bandwidth increases from 5 nm to 9 nm, the value of -Py^ sNstepOo^yp 4s decrease 



(2) 

to a value around 1. This decrease results mainly from the stronger increase of A\^l{u4s,3s) 
relative to ^sNstep@<^7p 4^ ('^4s,3s — 5) for 5 = and 5 ~ [see Fig. 5](b)]. Then, as the 
bandwidth further increases to around 15 nm, the value of -Pyp sNstepOwyp 4^ approach zero 



smce 



A 



(2+l)o7i— res 
7p,SNstep@a;7p,4s 



and 



A 



(2+l)near—res 
7p,SNstep@a;7p,4s 



0, resulting from A'^^^^^^^^^^ Ju^^ 



s,3s) 







and A 



(2) 



^SNstep@L^7p,4, (^4s,3s - 5) ^ for 5 ^ [scc abovc and Fig. sj^c)]. Actually, with 14.6- 
nm bandwidth the on-resonant component is completely zeroing out and to a very good 
approximation the corresponding shaped pulse is a three-photon dark pulse. 

Further increase of the bandwidth beyond 15 nm, leads to a continuous increase in 
the value of -Pyp sNstep@a;7p 4^ • "^^^ example, at 40-nm bandwidth it reaches a value of 
4.1. Following the above analysis and as seen in Fig. gdl)-(el) and Fig. ^a), for band- 
width values beyond 15 nm, the value of P^l^^step&ujrp 4^ actually determined only by the 



dominant near-resonant shaped-pulse term 



A, 



'^7p,4s 
(2+l)near— res 



A 



(2+l)on— res 
7p,TL 



ip p(2+l) 
, i-c., -1 7p^SNstep@aJ7p,4s ^ 



A 



(2+l)raear— res 
7p,SNstep@a;7p^4s 



7p,SNstep@ti)7p,4s 
2 



/ 



and the on-resonant TL term 
2 



A 



(2+l)on— res 
7p,TL 



. For large bandwidth 



Alu, E{uj7p^is) ~ E{uJo) = 1 and E^Ufp^^s — S) ^ E^luq) = 1 for small 6. Hence: (i) For the 

: kiAu, where ki is 



TL pulse: since a!^1{uj4s,3s) oc Alo, based on Eq. (3), 

a proportionality constant for the TL pulse; (ii) For the shaped pulse of tUgfgp'^*'''''' = uj7p^4^s- 
^7pS"tep@"J7p.4s = -pJZc^ lA^'^\uj4s,3s " S)d6 = k2aAuj + fcsfoAu; lu Au;, where k2a and ^26 
are constants for the specific single-step pulse shape. As a result, for large enough band- 
width values, an increase in the bandwidth Au leads to a continuous increase in the value 
-^7p SNstep@w7p 4s -^2a + -^26 1^ Acu. This is indeed seen in the graph shown in Fig. ij^a) 



A 



(2+l)o7i— res 
7p,TL 



(solid line). 
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^ {2+l)on—res 
7p,DBstep@(J7p,4s 



1 {2+l)on—res 



C. Shaped pulse with double-step phase pattern of "^[^t^step^ ~ '^7p,4s 

Last we consider the case of the shaped pulse with the phase pattern of a double vr step 
with a;[^^^g^*gp^ = ujp^^s (solid thin black lines in Fig. 5). With such a phase pattern, as with 
any of the double-step phase patterns, the phase associated with all the two-photon pathways 
of total two-photon transition frequency co'4s 3^ is zero, since $(0;) -|- $(co'4s^3s — tu) = for any 
u! (see Fig. [T]). So, fully constructive interferences are induced between all these two-photon 
pathways, as is the case with the transform-limited pulse. Hence, as seen in Fig. [5} with 
any spectral bandwidth, ^DBstep@w7p4s('^4s.3'') ~ ^tl('^4s,3s), which, as explained above, is 
the maximal possible value (here, at 5 = 0) for a given spectrum. This leads to a maximal 
on-resonant amplitude that is of equal magnitude to the TL case, i.e., 
, for any bandwidth Au. 
With the double-step phase pattern of t^ieft^^tep^ = ^7p,4s, the field -EDBstep@<^7p,4« (^7p,4s + 5) 
has different sign for ± \6\ when \6\ is small, which is in full correlation with the sign 
of 1/(5 for small ± This is similar to the case of the single-step phase pattern of 
'-^step'^*'^^'' = ^ipAs discussed above. Thus, also here, the relative sign between the ampli- 
tudes ^A^^^^l^^^^^^^^{6) contributed to the p-integration of 
determined only by ^DBstep(a^7p,4. (^4s,3s - S). 

The fact that ^DBstep@t<;7p4s('^4''.3«) (^•^•' ^'^^ S = 0) has the maximal possible value 

('2') 

also implies that for small 6 values ^DBstep@aj7p 4^ ("^^^.s* ~ ^) ^^^^ close to the 

maximal possible TL value, i.e., ^DBstep@<^7p,4, (^4s,3s - 5) ~ A\.l{ujis,3s - 5), with a cor- 
responding value that is a real positive one. As seen in Fig. [5| for any bandwidth, 

(2) 

^DBstep@a;7p 4s (^4«.3« ~ ^) '^^ Same sigu for small ± \6\, i.e., the interferences between 
negatively- and positively-detuned near-resonant pathways are constructive ones. Similar to 
the single-step phase pattern discussed above, the constructive nature of these interferences 
leads to the dominance of the corresponding near-resonant component ^7p^Bstep@^7p 4^ ^^^'^ 
the on-resonant component ^7p^B°tep@i!)7p 4^ bandwidth values. This is indeed the 

behavior of the near-resonant vs. on-resonant components shown in Fig. |4](b). The fact that 

(2) (2) 
^DBstepOtJTp 4s ('-^4«,3s ~ ~ ^TL ('^4s,3s — ^ ) , while the nature of the near-resonant interferences 

is constructive to the pulse with the double-step phase pattern and destructive for the TL 



< {2+l)near—res 
7p,DBstep@aj7p_4s 



2 



> 



A {2+l)near—res 
^7p,TL 



2 



pulse, imply that 

The overall bandwidth dependence of -^7^ DBstep@<^7p 4^ shown in Fig. ij^b) (solid line). 
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can then be understood as follows: Since 



A 



{2+l)on—res 
7p,DBstep@a;7p.4s 



A 



{2+l)on—res 
7p,TL 



and 



A 



{2+l)near—res 
7p,DBstep@aJ7p_4s 



> 



A 



(2+l)neaT—res 
7p,TL 



one obtains that P. 



(2+1) 



i.e., the corresponding three-photon 



width value. This is indeed what observed in Fig. |4 
absorption is always above the TL absorption. 

The asymptotic behavior of the graph of Pj^-Qstepmuiyp actually very similar 
to the single-step phase pattern case analyzed above, just with higher values due to the fact 



that A 



(2) 

DBstep@t<J7p,4s 



(i^4s,3s) is kept on its maximal value for any bandwidth. Also here, for 



large bandwidth Atu, following the above analysis for the present double-phase-step pulse 



and for the TL pulse, P^p_DBstep 



A. 



(2+l)near—res 
7p,DBstep@cJ7p_4s 



/ 



A 



(2+l)o7i— res 
7p,TL 



. Also, E(u;7p,4s) 



(2') 

E{u!o) = 1 and E{u'rp^4^s — S) ~ E{ujo) = 1 for small 6. So: (i) For the TL pulse: Axl(^4s,3s) ^ 

kiAuj, where ki is the TL proportionality constant, (ii) For 



Auj and thus 



A, 



(2+l)on— res 
7p,TL 



,11 1 1 r (DBstep) < (2+l)near— res 

tne snaped pulse oi ^i;^ 

-steD ^7p,4s- ^7p,DBstep^ 



^left-step ~ '^Yp,4s- ^ ■^7p,DBstep@<^7p,4s ~ 5" 

ksaAuj + k^bAuj \n Auj , where k^a and fcsf, are constants for the specific double-phase-step 
pulse shape. So, -^7^ DBstep@<^7p 4^ behaves asymptotically as K^a + -ft'sftln Acu. The constants 
K^a and K^^ are, respectively, larger than the and constants of the single-step 

'^7p,4s5 due to the larger near-resonant component induced by 
the double-step phase pattern. Hence, for any given bandwidth, the resonance-mediated 
(2+1) three-photon absorption induced by the shaped pulse with double-step phase pattern 
of i-<^w+^ff™^ = '^7p,4s is higher than the absorption induced by the shaped pulse with single- 



phase pattern of i^itep^*^''^ 



left-step 

step phase pattern of cUgf^p'^*'^^'' 



l^7p,4s- 



V. CONCLUSIONS 



In conclusion, we have studied in detail femtosecond phase control of atomic resonance- 
mediated (2+1) three-photon absorption, focusing on its dependence on the spectral band- 
width of the excitation pulse. The total absorption amplitude has contributions from a 
group of interfering three-photon pathways that are on resonant with the intermediate state 
and a group of interfering three-photon pathways that are near resonant with it. For a given 
pulse shape, the ratio between the amphtude contributions from these two groups depends 
on the spectral excitation bandwidth. The nature of the three-photon absorption that is 
induced by a given shaped pulse is analyzed in a systematic and physically-intuitive way 
by decomposing it into a two-photon absorption, which is induced by all the possible two- 
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photon pathways leading to various detunings from the intermediate state, followed by the 
absorption of a proper third photon. Prominent among the results is our finding that with 
simple proper pulse shaping an increase in the excitation bandwidth leads to a corresponding 
increase in the enhancement of the three-photon absorption over the absorption induced by 
the (unshaped) transform-limited pulse. For example, here, 40-nm bandwidth leads to an 
order-of-magnitude enhancement over the transform-limited absorption. The present work 
serves as a basis for future extensions to molecular systems and we expect it to be significant 
for nonlinear spectroscopy and microscopy. 
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FIG. 1: The resonance-mediated (2+1) three-photon excitation scheme (not to scale) from j^f) 
to I/) via \r). The corresponding Na states are, respectively, 3s, 4s, and 7p. Several sets of 
pathways that are on resonance {6 = 0) and near resonance {6 ^ 0) with |r) = 4s are shown; 6 is 
the corresponding detuning. Also shown schematically are examples of the two types of spectral 
phase patterns used in the study, $(SNstep)(^) ^(DBstep)^^) The $(SNstep)(^) includes the 
phase patterns with a single tt step, each is characterized by the vr step position i-i-'g^gp ■ The 
^(DBstcp)^^-j set includes the phase patterns composed of two vr steps, with the two steps located at 
equal distance along the frequency axis from both sides of half the two-photon transition frequency 
ojrgf^- Each such double-step pattern is characterized by the position of the left step (^i^t^g^tep ^ • 
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FIG. 2: Experimental (circles and squares) and numerical-theoretical (solid lines) results for the 
resonance-mediated (2-1-1) three-photon absorption in Na. Panels (al) and (a2) correspond to 
the full pulse spectrum [inset of (al)], while panels (bl) and (b2) correspond to the blocked pulse 
spectrum [inset of (bl)]. Panels (al) and (bl) (left-column panels) present results for the single-step 
set of spectral phase patterns $(S'^^*^P)(a;). The corresponding traces show the TL-normalized final 
7p population i^y^^^'' (see text) as a function of the single vr-step position Wgfgp^*'^'''' . Panels (a2) 
and (b2) (right-column panels) present results for the double-step set of spectral phase patterns 
^(DBstep)^^-j The corresponding traces show the TL-normalized final 7p population i-yp^^"* (see 
text) as a function of the left 7r-step position oj^^^^ . 
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FIG. 3: Theoretical results (solid lines) of the resonance-mediated (2+1) three-photon absorp- 
tion in Na calculated with excitation pulse of different spectral intensity bandwidth Aco. In all 
the cases the pulse spectrum is a Gaussian centered at 780 nm (12821 cm~^). Each row in the 
figure corresponds to a different bandwidth. The left-column panels correspond to the single-step 
spectral phase patterns $(S^^*^P)(a;). They show the TL-normalized final 7p population Pjp^^^ 
(solid lines) as a function of the single 7r-step position cj^f^p'^^'^^^ Additionally, they also show the 



corresponding on-rcsonant component 

2 

'T(2+l)near—res 
^7p 



2 

,{2+l)on—res 
^7p 



(dotted lines) and ncar-rcsonant component 



(dashed lines). The right-column panels correspond to the double-step spectral 
phase patterns #(°^^*^P)(a;). They show the TL-normalized final 7p population P^p^^^ as a function 
of the left TT-step position i^i^^^l^p^ ■ 
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FIG. 4: Numerical-theoretical results for the dependence of the resonance-mediated (2+1) three- 
photon absorption in Na on the spectral intensity bandwidth of the excitation pulse having (a) 
single-step phase pattern of Wgf^p****''''* = a;7p,4s, and (b) doublc-stcp phase pattern of ^[^f^^^^^ = 



^7p,4s- Shown are pj^'*'^'* (solid lines), 
lines) as a function of the spectral bandwidth Aa;. 
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FIG. 5: Numerical-theoretical results for the values of A^'^\Q = a;4s,3s — 5) [see Eq. ^] as a 
function of the detuning 6 for different spectral intensity bandwidth Acu of the excitation pulse. 
In all the cases the pulse spectrum is a Gaussian centered at 780 nm (12821 cm^^). The different 
panels correspond to the different panels in Fig. [3j The detuning values are given as normalized 
values 6/Auj (see text). Each panel shows results for the transform- limited (TL) pulse (dashed 
black line), for the shaped pulse with the single-step phase pattern of co'sfep***'^''^ = ^^^7^,45 (solid gray 
line) , and for the shaped pulse with the double-step phase pattern of "^i^t^step ^ ~ '^7p,4s (solid black 
line) . 
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